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ABSTRACT 

We present optical morphologies obtained from deep Hubble Space Telescope {HST) and ground-based 
images for galaxies selected from the first sub-millimeter (sub-mm) survey of the distant Universe. Our 
sample comprises galaxies detected in deep 850-/im continuum maps of seven massive clusters, obtained 
using SCUBA, the new bolometer camera on the JCMT. The survey covers a total area of 0.01 degree^ 
to Icr noise levels of about 2mJybeam~^. We detect a total of 25 sources at 850 /xm, of which 17 
and 10 are brighter than the respective 50% and 80% completeness limits. Optical counterparts are 
identified for 14 of the 16 sources in the f^Q% sample and for 9 of the 10 sources in the fsQ% sample 
that lie within our optical fields. The morphologies of those galaxies for which we have HST imaging 
fall into three broad categories: faint disturbed galaxies and interactions; faint galaxies too compact to 
classify reliably; and dusty, star-forming galaxies at intermediate redshifts. The disturbed and interacting 
galaxies constitute the largest class, which suggests that interactions remain an important mechanism 
for triggering star formation and the formation of ultraluminous galaxies in the distant Universe. The 
faint, compact galaxies may represent a later evolutionary stage in these mergers, or more centrally- 
concentrated starbursts. It is likely that some of these will host active galactic nuclei. Analysis of the 
colors of our sample allow us to estimate a crude redshift distribution: > 75% have z < 5.5 whilst 
> 50% lie at z < 4.5, suggesting that the luminous sub-mm population is coeval with the more modestly 
star-forming galaxies selected by UV/optical surveys of the distant Universe. This imposes important 
constraints on models of galaxy formation and evolution. 

Subject headings: cosmology: observations — cosmology: early universe — galaxies: evolution — 
galaxies: formation — galaxies: morphology 



1. INTRODUCTION 

Recent technological advances in ground-based sub-mm 
technology (Holland et al. 1998) are set to revolutionise 
our understanding of the star-formation history of galax- 
ies (Blain et al. 1998a, BSIK). Deep maps of the sub- 
mm sky (Small, Ivison & Blain 1997, SIB) have now re- 
solved the population of galaxies responsible for the far- 
infrared background (Puget et al. 1996; Fixsen et al. 1998; 
Schlegel, Finkbeiner & Davis 1998) and offer a precise de- 
termination of the star- formation rate (SFR) in distant, 
dust-obscured systems (e.g. Dey et al. 1998). This popu- 
lation is responsible for a substantial fraction of the total 
SFR in the early Universe (BSIK) and holds crucial infor- 
mation about how galaxies formed and evolved (Blain & 
Longair 1993). 

Our sub-mm survey, the first results of which were re- 
ported by SIB, exploits massive, foreground clusters to en- 
hance the sensitivity of the maps to distant star-forming 
galaxies (Blain 1997). The survey has now been extended 
to cover seven clusters (Small et al. 1998b, S98) with suffi- 
cient sensitivity to detect a luminous infrared galaxy such 
as Arp220 out to z ~ 10. Our sub-mm galaxy counts 
greatly exceed those expected in a non-evolving model 
based on the local IRAS 60-fj,m luminosity function (Saun- 
ders et al. 1990), or in a model that can explain the 

^Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy Inc., under NASA contract NAS 5-26555. 

3PPARC Advanced Fellow. 
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optically-selected Lyman-dropout (Madau et al. 1996) and 
Lyman-emission samples (Hu, Cowie & McMahon 1998). 

As demonstrated by BSIK, before we can understand 
the star-formation history of the Universe, we need one 
remaining piece of the puzzle: the redshift distribution, 
N{z)^ of a faint sub-mm-selected galaxy sample. The iden- 
tification of the complete S98 catalog is underway (Ivi- 
son et al. 1998a), this benefits from archival 775' T" images 
for most of the survey area as well as the availability of 
10-/xJy beam~^ radio maps for all the fields (Ivison et al. 
1998b). In the meantime, useful limits on the sub-mm 
N{z) can be determined from broad-band imaging (Stei- 
del et al. 1996). Even a basic analysis along these lines can 
challenge some currently viable models of galaxy evolution 
(BSIK). 

In this paper we present limits on the redshift distri- 
bution of a sample of sub-mm-selected galaxies derived 
from deep optical imaging. We also take advantage of 
archival HST images to investigate the optical morpholo- 
gies for a sub-sample of our catalog on angular scales down 
to O.larcsec (< 1/i^Q^kpc at the expected redshifts). This 
detailed morphological information can provide a valuable 
insight into the physical mechanisms responsible for the 
high SFR of these distant galaxies and thereby providing 
a fuller understanding of the early stages of galaxy forma- 
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Table 1 

Log of the SCUBA observations 
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of these images is documented in those papers. While 
lacking the spatial resolution and depth of the WFPC2 
images, these data have a wider field of view and a more 
comprehensive wavelength coverage. 

Total magnitudes for all the sub-mm-selected galaxies 
covered by our optical images were estimated using SEx- 
tractor's BEST_MAG in the reddest passband available 
for each cluster (either Ror I), using ffST frames if avail- 
able. Galaxy colors were measured within 1-arcsec diame- 
ter apertures on the WFPC2 frames and 3-arcsec apertures 
on the ground-based images, after matching the effective 
seeing on the different frames. If required, the photometry 
was converted to the /-band assuming (R — I) ~0.5±0.2, 
which is typical of both the faint field population (Small 
et al. 1995) and the available photometry for our sub- 
mm counterparts. For 010939+47, the {g — i) colors from 
Dressier & Gunn (1992) were transformed to (V — I) (Kent 
1985). No reddening corrections were applied. 

2.3. Optical Identifications 

The absolute positional accuracy of the SCUBA maps 
is better than 3 arcsec, as estimated from repeated point- 
ing checks performed during the observations. The posi- 
tions of detected sources were determined by fitting a cir- 
cular gaussian to a 30-arcsec diameter region (two beam 
widths) around the peak flux position on 4-arcsec-sampled 
maps. The random errors in the positions of individual 
objects depend on the signal-to-noise ratio of the detec- 
tions, but should be no worse than 4 arcsec rms for the 
faintest sources discussed here. The combined positional 
uncertainty of our faintest sub-mm sources is therefore 
^ 6 arcsec (c.f. Fig. 1, Plate. 1). Counterparts were sought 
within this radius of the nominal sub-mm position. 

Using the deepest and/or reddest passband available for 
each cluster, we identified all counterparts down to the 
magnitude limits of our optical data for the 25 sub-mm 
sources from the full catalog. These limits are typically 
/ ~ 23.5 and 26.0 for the ground-based and i/S'T exposures 
respectively. Ten sources in the /8o% sample lie within our 
optical images; 7 of these are on the HST frames and the 
remaining 3 on ground-based images. The equivalent num- 
bers for the 16 source in the f^Q% sample are 11 and 5. In 
two maps a bright sub-mm source has a clear counterpart, 
and so their coordinates were corrected to align the sources 
exactly. 

Due to the nature of sub-mm-selected galaxies, the opti- 
cal coimterparts are likely to be faint, moreover the cluster 
fields analysed here are crowded. Thus in some instances 
there is a non-zero likelihood of an unrelated galaxy falling 

*The JCMT is operated by the Observatories on behalf of the UK Particle Physics and Astronomy Research Council, the Netherlands 
Organization for Scientific Research and the Canadian National Research Council. 



2. OBSERVATIONS, REDUCTION AND ANALYSIS 

2.1. Suh-millimeter Observations 

Our 850-/im maps were constructed using the long- 
wavelength array of the Sub-millimeter Oommon-User 

Bolometer Array (SCUBA, Robson et al. 1998) on the 
James Clerk Maxwell Telescope (JCMT)^. The final on- 
source integration times and field centers are listed in Ta- 
ble 1. Each field covers an area of 5.2 arcmin^, giving a to- 
tal survey area of 0.01 degree^. The typical amplification 
for background sources detected in our fields is expected 
to be 1.5to.5) ^'^'i have in efi'ect surveyed an area of 
about 25 arcmin^ in the source plane to an equivalent Icr 
sensitivity of 1.3 mJy. 

Source catalogs were constructed using the SExtractor 
package (Bertin & Arnouts 1996) and Monte Carlo incom- 
pleteness simulations undertaken as described in SIB, for 
more details see Blain et al. (1998b) and S98. The derived 
sensitivity limits, 850-/xm fluxes {fs5a) for 50% and 80% 
completeness of /85o(50%) and /85o(80%), and the number 
of detected sources in each field are listed in Table 1. We 
identify only two sources by applying the same selection 
criteria to the negative fluctuations in our maps, both have 
/850 < ./85o(50%), indicating a false detection rate in our 
full catalog of ^ 10%. In the remainder of our analysis we 
will focus on the two sub-samples from our full catalog of 
25 sources: 10 sources with /gso > /85o(80%) which we will 
denote /8o% and a slightly more liberal sub-sample (/5o%) 
defined by fs5o > /85o(50%) and containing 17 sources. 
The /8o% and /5o% limits are roughly equivalent to 4- and 
3- (7 of the sky noise. 

2.2. Archival Optical Observations 

In Table 2 we hst the HST WFPC2 and the wider- 
field, ground-based observations used in our analysis. 
The WFPC2 observations are mainly from programs with 
which we are associated; the remainder were retrieved from 
the ST-ECF archive. All were reduced as described by 
Small et al. (1997). For one cluster, C10939-h47, we have 
used the deconvolved pre-refurbishment WF/PC-1 images 
of Dressier et al. (1994) to determine the morphologies for 
SCUBA sources lying outside the WFPC2 field of view. 
The ground-based observations come from a variety of 
sources, predominantly from programs on the 3.6-m CFHT 
(Kneib et al. 1993) and the 5.1-m Hale (Dressier & Gunn 
1992; Small et al. 1998a). The reduction and calibration 
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Table 2 

Log of archival optical observations 



Cluster 


HST/WFPC2 
(Filter [Exposure timc/ks]) 


Reference 


Ground-based 


Reference 


CI 0024+16 
A370 

MS 0440+02 

CI 0939+47 

A1835 

A2390 

CI 2244-02 


F450W [23.4] F814W [13.2] 
F336W [27.8] F675W [5.6] 
F702W[18.4] 
F702W[21.0] F702W[22.0]° 

F555W [8.4] F814W [10.5] 
F555W [8.4] F814W [10.5] 


Smail et al. (1997) 

Archive 

Archive 

Smail et al. (1997) 

Kneib et al. (1998) 
Archive 


/ 

UBVRI 

gri 
UBI 
UBI 


Mellier (priv. comm.) 
Kneib et al. (1993) 

Dressier & Gunn (1992) 
Smail et al. (1998) 
Smail et al. (1998) 



a) WF/PC-1 observations, reduced and deconvolved as described in Dressier et al. (1994) 



within the sub-mm error box. As such our identifications 
should be treated as preUminary until they have been con- 
firmed from deep radio maps, in particular it is possible 
that in a small number of cases we have misidentified a 
sub-mm source due to the faintness of the true counterpart 
combined with the presence of a nearby, brighter galaxy. 
Nevertheless, as we show below the majority of the iden- 
tified counterparts are likely to be real and so the broad 
characteristics of luminous sub-mm galaxies can be ascer- 
tained from our current catalog. 

The likelihood of a false identification for a given sub- 
mm source will be a function of both the apparent mag- 
nitude of the identified counterpart and its distance from 
the nominal position of the sub-mm source. We have used 
this information to estimate the probability of a real as- 
sociation on an objcct-by-objcct basis. Using a Monte 
Carlo simulation, based upon the observed number counts 
of galaxies within each of our optical images, we estimate 
the probability of a galaxy with the observed apparent 
magnitude (or brighter) falling at random within a circle 
defined by the optical and sub-mm positions. The reliabil- 
ity of an identification have been assigned to one of four 
classes based upon the probability of a chance coincidence 
P: P < 0.01, definite; P < 0.10, likely; P < 0.25, possible 
and P > 0.25, blank field. 

3. RESULTS AND DISCUSSION 

The aim of this work is to obtain a first view of the 
optical properties of the galaxies which contribute a sub- 
stantial fraction of the sub-mm emission in the Universe. 
To this end we retain the /5o% sample although it is not 
statistically complete, instead we use it to provide support 
for any trends seen in the smaller /8o% sample. In the /8o% 
sample of 10 sources within the optical fields, there are 4 
counterparts with P < 0.0, 9 with P < 0.25 and 1 blank 
field. The rms offset between sub-mm sources and the 
9 counterparts is 2.3 ± 0.5 arcsec, in reasonable agreement 
with the nominal positional accuracy of the SCUBA maps. 
The 16 sources with /gso > /85o(50%) are associated with 
6 definite identifications, 8 likely/possible identifications 
and 2 blank fields. 

In the following discussion we use all counterparts with 
P < 0.25, noting that none of our conclusions arc sub- 
stantially altered if we restrict our sample to those coun- 
terparts with P < 0.10. 

Multicolor photometry is available for 9 of the 10 sources 
in the fgo% sample (12/16 in the /8o%)- One of these 9 
fields is blank and in the remaining 8 we have 6 B and 
2 V-band detections. For the more liberal /5o% sample 
we have 8 B-band and 4 V-band detections from a total of 



12 identified sources with multicolor photometry. In all in- 
stances where we have both B and F-band data, the coun- 
terpart galaxies are detected in both bands. The galaxies 
plotted in Fig. 2 have a broad apparent-magnitude distri- 
bution and span the whole range of colors seen in the faint 
field population. We find no strong correlation between 
sub-mm fiux and /-band magnitude, as expected due to 
the wider selection function of a sub-mm survey compared 
with that of optical surveys. 

The detection of a galaxy in one of the bluer passbands 
can be used to impose an upper limit to its rcdshift, be- 
cause of the strong absorption by the foreground Lyman-a 
forest at rest-frame wavelengths shortward of 912A. Thus 
the number of detections in the B- and y-band images 
requires that at least 8/10 or > 80% of the galaxies in 
our /8o% sample have z < 5.5, the rcdshift at which the 
Lyman-limit moves through the ^-band. At least half of 
this sample lies at z ^ 4.5 based on the B-band data alone. 
The equivalent limits for the /5o% sample are ^ 75% at 
z < 5.5 and > 50% with z < 4.5. Thus these luminous 
sub-mm galaxies appear to span a similar rcdshift range 
to that covered by the more modestly star-forming galax- 
ies selected in UV/ optical surveys of the distant Universe 
(Madau et al. 1996). 
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Fig. 2. The {V - color-magnitude diagram for the 18 SCUBA 

sources for which wc have muKi-color information, the approximate Um- 
its for the 2 non-detections are marked. The filled symbols and sohd 
limits show those sources in the f80% sample, the open points and dot- 
ted limits are the remaining /5o% sources. The square symbols represent 
{B — I) measurements converted to {V — I) assuming {B — V) — 0.5. The 
apparent magnitudes of the background sources have been corrected for 
a mean amplification of 1.5. The dashed line shows the effective com- 
pleteness limit in our optical images. The solid line indicates the median 
{V — I) color as a function of I magnitude for the faint field population, 
and the histogram gives the distribution of {V — I) colors for a sample 
of field galaxies down to / = 25 (Smail et al. 1995). 

The range of evolutionary models that can account for 
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all the sub-mm and far-infrared data currently available 
(BSIK) can already be constrained by these crude hmits. 
We can reject any model which predicts that more than 
~ 20% of galaxies selected at an 850- /im flux limit of 3 mJy 
lie beyond z ~ 5.5. 

Turning to the morphological composition of our cat- 
alog, we have used high-quality HST images to provide 
detailed morphologies for the 11 sub-mm-selected galaxies 
within these fields (Fig. 1). For the remaining 5 objects we 
use ground-based images, with typically sub-arcsec resolu- 
tion, to make cruder classifications. Across the whole cat- 
alog we find that the morphologies of the sub-mm galaxies 
can be split roughly into 3 classes: 1) faint, disturbed 
galaxies, some of which are obviously undergoing mergers 
and interactions; 2) faint, compact galaxies, too small to 
classify reliably; and 3) bright, dusty star-forming galaxies. 
We define a disturbed galaxy on the basis of the presence 
of one or more of the following; multiple components, close 
companions with similar sizes and colors, clearly distorted 
morphologies or tidal extensions. Some of these signatures 
are only visible on scales of ^ 1 arcsec in the RST images, 
and so there is an obvious bias against identifying such sys- 
tems in our ground-based data. The morphologies of the 
9 identified counterparts in the /8o% sample fall into the 
three classes in the ratio 5(2):1(1):3(0), numbers in paren- 
theses give the values from the ground-based observations. 
The corresponding ratios for the 14 optical counterparts 
in the /850 > /85o(50%) sample are 6(3):4(1):4(1). 

The bright galaxies, shown in the bottom row of Fig. 1, 
include two central galaxies from the clusters used in our 
survey. These clusters both contain strong cooling flows 
and their central galaxies show signatures of strong star- 
formation activity and some non-thermal emission. The 
properties of these two galaxies are investigated in more 
depth in a forthcoming paper. As confirmed cluster mem- 
bers these two galaxies are removed from the analysis here. 
The remaining 2 bright galaxies are both dusty Sab~Sbc 
spirals, lying at intermediate redshifts, within or in front 
of the clusters. The flux density limit of our 850-/im sur- 
vey means that we would expect to detect this type of 
galaxy, with a SFR > 20-50 Mgyr"!, at z < 0.2-0.5. We 
note that none of these bright galaxies was included in the 
analysis of SIB. 

The optically faint sources are most interesting, as 
they typically lie at high redshifts and have high SFRs, 
> lOOM0yr~^. The number of these galaxies that are 
clearly disturbed/interacting systems is roughly equal to 
the number that are either too faint or too compact to be 
classified reliably. This is a very robust result - at least 
50% of our faint, distant sub-mm-selected sources are iden- 
tified with distorted or interacting galaxies. This interac- 
tion rate is higher than the 30% seen in the general field 
population at a similar magnitude limit (van den Bergh et 
al. 1996). The median separation of the observed compo- 
nents in the interacting systems is about 10/i^q^ kpc. 50% 



is thus a firm lower limit to the fraction of interactions 
as some of the more compact sources, with typical half- 
light radii of 0.2-0.4 arcsec, or 2-4ft,^Q^ kpc at z > 1, may 
be either unresolved or more evolved mergers. The exact 
nature of the compact galaxies will require further investi- 
gation, in particular we are obtaining high-resolution near- 
infrared images to provide a more representative view of 
these galaxies. 

Strong interactions and mergers are mechanisms respon- 
sible for the high SFRs of ultraluminous infrared galaxies 
(ULIRG) in the local Universe (Sanders & Mirabel 1996). 
The fraction of interacting ULIRGs found locally is compa- 
rable to or greater than our lower limit, but the mean sep- 
aration of the components is close to the resolution limit of 
even our WFPC2 images. Hence we suggest that the rate 
of interaction found in our distant, sub-mm-selected galax- 
ies is probably comparable to that seen in local ULIRGs. 
On the basis of our morphological survey we can state that 
one mechanism, dynamical interaction, is responsible for 
triggering at least half of the sub-mm flux density, and by 
implication for at least half of the star-formation activity 
in the most distant strongly star-forming galaxies in our 
sample. These observations provide the clearest support 
for the pivotal role of mergers and interactions in the early 
evolution of luminous star- forming galaxies. 

4. CONCLUSIONS 

• We present optical identifications of galaxies selected in a 
deep sub-mm survey. Roughly 70% of the sources selected 
at 850 /ini, for which optical imaging data is available, have 
definite/likely counterparts. A further 20% of the sub-mm 
sources have probable optical counterparts, only 20% of 
the sub-mm sources have no identifiable counterpart to a 
typical apparent magnitude of / ~ 25-26. 

• 75% of the optical counterparts for which we have B 
or y-band images are detected in these passbands, sug- 
gesting that at least three-quarters of the 850-^m-selected 
galaxies brighter than 3mJy lie at z ^ 5.5. 

• The morphologies of the fainter (and hence probably 
more distant) optical counterparts fall roughly equally into 
two classes: disturbed/interacting and compact objects. 
The prevalence of interactions in this survey underlines 
the central role of this mechanism in triggering star for- 
mation in the most luminous galaxies at high redshift. 
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Fig. 1. 10 X lOarcscc images of the 16 sub-mm sources in our /5o% sample whieh fall within our optieal imaging of the seven clusters. These are 
ordered from the upper-left on the basis of their morphologies: 6 disturbed/interacting, 4 compact/featureless (including a strongly-distorted arclet), 
2 blank fields and then 4 bright, dusty galaxies. A '?' next to the sub-mni source name indicates a blank field, a marks those sub-mm sources for 
which /sso > /s5o(80%)- The images are centred on the most likely optical candidate. The centroids of the sub-mm sources are indicated by crosses. 
Note that these images span a range in exposure times and resolutions (Table 2). The panels correspond to ^ 80h^Q kpc at 2 > 1. 



